Epitaxial growth of binary alloy nanostructures 



S. Heinrichs, 1 W. Dieterich, 1 and P. Maass 2 
1 Fachbereich Physik, Universitdt Konstanz D-78457 Konstanz, Germany and 
2 Technische Universitdt Ilmenau - Institut fur Physik, D-98684 Ilmenau, Germany 

(Dated: July 11, 2006) 

Stochastic growth of binary alloys on a weakly interacting substrate is studied by kinetic Monte 
Carlo simulation. The underlying lattice model relates to fee alloys, and the kinetics are based 
on deposition, atomic migration with bond-breaking processes and exchange processes mediated by 
nearest neighbor hopping steps. We investigate the interrelation between surface processes and the 
emerging nonequilibrium structure at and below the growing surface under conditions where atoms 
in the bulk can be regarded as immobile. The parameters of the model are adapted to CoPt3 alloys. 
Growing nanoclusters exhibit an anisotropic short range order, primarily caused by Pt segregation 
at the surface. The overall structural anisotropy depends on both Pt surface segregation and cluster 
shape, and can explain the perpendicular magnetic anisotropy (PMA) recently measured in CoPt3 
nanoclusters on a van der Waals substrate. The onset of LI2 ordering in the cluster is induced by 
surface processes. The same kinetic model is applied also to continuous thin films, which in addition 
can exhibit a small bulk contribution to PMA. 

PACS numbers: 81.15.Aa,68.55.Ac 



I. INTRODUCTION 

Molecular beam epitaxy (MBE) has become an impor- 
tant tool to prepare ultrathin films and nanostructures in 
nonequilibrium, frozen-in states that show novel proper- 
ties distinctly different from the equilibrium bulk phases. 
Theoretically, an important question is how to relate the 
structural characteristics of the growing material to the 
incoming flux, substrate temperature, adatom diffusion 
coefficients and adatom interactions. Successive stages 
of growth, nucleation and island formation in the sub- 
monolayer regimeji second layer nucleationj 2 ^ growth of 
three-dimensional clusters and emerging film morpholo- 
gies, are fairly well understood for the one-component 
case, i.e. for deposition of only one species of atoms4£ By 
contrast, nonequilibrium alloy nanostructures produced 
by codeposition of two or more atomic species can dis- 
play a variety of additional new phenomena that are only 
partly understood up to now. From the viewpoint of the- 
oretical modelling, this is true in particular for metallic 
nanoalloys, which recently became of great interest from 
the experimental and technological viewpoint & 

In many MBE experiments the substrate temperature 
is chosen sufficiently low so that atomic configurations in 
the interior of a growing cluster or film can be considered 
as frozen. Structural relaxation then takes place only 
among lower-coordinated atoms within the growth zone. 
Specific to binary systems are the following questions: 

i) Suppose that one atomic species tends to segregate 
at the surface. What is the emerging bulk struc- 
ture upon further deposition burying a segregated 
surface layer? 

ii) How large is a possible anisotropy in the frozen 
short-range order, i.e. a difference in the alloy struc- 
ture between the lateral and the perpendicular 
(growth) direction? 



iii) To what extent can diffusion processes limited to 
the near-surface region give rise to long-range order 
in the bulk of an ordering alloy at temperatures be- 
low the equilibrium order-disorder transition point? 

iv) In the case of magnetic alloys, what are the 
magnetic properties associated with the non- 
equilibrium short or long-range order? This ques- 
tion is central in the ongoing search for materi- 
als that display a stable perpendicular magnetic 
anisotropy (PMA), which is useful for the develop- 
ment of high-density magnetic storage media. Un- 
der this viewpoint fcc-alloy systems of Co-Ptj^ Fe- 
PliS or Fe-CoiS draw great attention in current in- 
vestigations. 

Here we report on answers related to these problems, 
based on kinetic Monte Carlo (KMC) simulations of a 
statistical model for binary fee alloys of composition AB3, 
with emphasis on CoPt3 alloys. 

MBE-grown nanoclusters of CoPt3 on a van der 
Waals WSe2 substrate were found recently to exhibit 
PMA at room temperature and below the onset of LI2- 
orderingi The origin of PMA is thought to lie in a 
structural anisotropy induced by the growth process, 
strong enough to overcome the dipolar form anisotropy 
favoring in-plane magnetization. The model we use is 
based on effective chemical interaction and kinetic pa- 
rameters extracted, respectively, from known equilib- 
rium properties of CoPtgii*^ and from separate diffu- 
sion measurements^ The dependence of the magnetic 
anisotropy on the local structure is represented by a 
sum over bond contributions, with nearest neighbor Co- 
Pt bonds as the dominant part. The associated bond 
anisotropy parameter is deduced from magnetic measure- 
ments on Co-Pt multilayer systems^ Within this model 
it was shown that Pt surface segregation and cluster 
shape effects can explain the occurence of PMA in some 
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temperature window, followed by the evolution of LI 2- 
order at higher temperatures^ In this paper we present 
further details on cluster structures and extend our pre- 
vious investigations to homogeneous films. 

In section UTI we introduce our kinetic model, based on 
random deposition and subsequent hopping moves of A 
and B atoms to neighboring sites of an fee lattice. Hop- 
ping rates are derived from the energetics in the specific 
environment of each atom. To our knowledge, most prior 
KMC simulations of three-dimensional growth of binary 
systems were in the framework of solid on solid (SOS) 
models^ which are not well suitable to describe cluster- 
shapes and facetting in a realistic way. In very recent 
work the growth of CoPt3 films was simulated, using a 
more realistic tight binding, second moment approxima- 
tion to the interatomic interaction potential^ Segrega- 
tion of Co-atoms to step edges on the growing surface 
was predicted, which leads to in-plane Co-clustering. 1 ^ 
By this the authors were able to explain PMA in dis- 
ordered films, known from previous measurements^ to 
occur at elevated temperatures (T > 500K). 

Our model, based on effective nearest neighbor cou- 
plings, allows for fully three-dimensional kinetics, a fea- 
ture which is important in simulations of cluster growth 
and LI 2 ordering. It yields cluster facets in agreement 
with those observed in experiments and consistently ac- 
counts for mass transport along the cluster side and 
top facets. Both height and lateral growth largely re- 
lies on interlayer diffusion that starts from adatoms on 
the substrate, reaching higher layers along side facets. 
The present model automatically entails step edge bar- 
riers, dissociation from the surface and the possibility of 
vacancy diffusion in the bulk. It is less comprehensive 
than recent simulations of island formation in compound 
semiconductors, which are based on extensive electronic 
density functional calculations for energy barriers of the 
processes involved. These calculations, however, were so 
far limited to the submonolayer regime^ Section III Bl 
also describes a phenomenological relationship between 
the local contributions to the magnetic anisotropy and 
atomic short range order within a bond picture. 

The anisotropic short range order and associated mag- 
netic properties of clusters up to 4000 atoms are pre- 
sented in section ITTTI and are discussed in detail in terms 
of surface segregation and cluster shape. The same algo- 
rithm is applied in section llVl to thin films. In addition 
to surface -induced PMA a bulk contribution is identified 
which, however, is too small to account for the measured 
PMA in thick films. Further conclusions are drawn in 
section Ivl 



II. SIMULATION MODEL 

A. Interactions and Kinetics 

Our model starts from a reference fee lattice with 
nearest neighbor distance a in a cubic simulation box 



L x L x L z where sites can be occupied by an atom of 
type A (Co), B (Pt), or by a vacancy V. In contrast to 
models with a small concentration of vacancies and fixed 
shape of the sample, most of the simulation box con- 
sidered here consists of vacancies representing the free 
space in the vicinity of the cluster. This allows for an 
unconstrained evolution of the cluster morphology. Ef- 
fective interactions between A and B atoms are restricted 
to nearest neighbor pairs. This simplified description al- 
ready captures essential statistical properties of ordering 
fee alloys, including the phase diagram which displays 
LI2 and Llo ordered structures in the vicinity of the AB3 
and AB-composition, respectively>i2i2S 

Bond energies between the different atomic species are 
denoted by VaAjVab and Vbb- The linear combina- 
tions I = j(Vaa + Vbb - 2Vab), h = Vbb - Vaa which 
acts as a surface field, and Vb = (Vab + Vbb) /2 control 
the bulk order-disorder transition temperature, the de- 
gree of surface segregation of B-atoms, and the average 
bond energy in the LI2 ordered state. Parameter val- 
ues are adjusted to reproduce equilibrium properties of 
CoPt3: (i) The transition from the disordered to the LI2- 
structure occurs at T ~ 958 K, 21 which is related to I by 
k B T ~ 1.837S Setting 7=1, our energy unit becomes 
fc s T /1.83 = 45meV, corresponding to 523 K. (ii) The 
observed Pt surface segregation of nearly 100%jii caused 
by the larger size of Pt relative to Co, is compatible with 
h > 4; here we choose h = 4 in most of our calculations, 
(iii) At temperatures of interest, the mobility of atoms is 
effectively restricted to the film or cluster surface so that 
their typical coordination is between 3 (for an adatom 
on top of a terrace) and 7 (for an atom attached to a 
step edge). Parameters for a variety of corresponding 
processes were calculated for Pt within electronic den- 
sity functional theory (DFT)^2, Different processes were 
classified according to the number of broken bonds, and 
the average bond energy yielded Vb ~ —5, which we use 
here. Clearly, these bond energies are to be understood 
as effective energies since real interactions can extend to 
several neighbors. 24 

Experiments for nanoclusters of CoPt3 were performed 
on a WSe2 (0001) substrate surface^ which is of the van 
der Waals type. The interaction with the substrate is 
modelled by a weak attractive potential, represented by 
an additional energy V s for atoms in the first layer. In 
the cluster simulations we choose V s = —5. Compared 
to an fee Pt (111) surface with three bonds of typical 
strength Vb, this amounts to about 1/3 of the energy of 
a single Pt-Pt bond. 

The total flux F of incoming atoms is taken as F — 3.5 
monolaycrs/s. Additional parameters needed to describe 
the kinetics of the system are the transition state energy 
Ut and the attempt frequency v for atomic moves. These 
values are known from diffusion experiments 1 ^ of Pt on 
Pt (111) surfaces as v = 8.3 • 10 11 s" 1 and U t ~ 5. The 
resulting diffusion coefficient is D/a 2 = {v / 4)e~ Ut l kBT 
for moves without change in energy. 

The rate for an elementary hopping process with a final 
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energy Ef and initial energy E{ is chosen to be 

Wif = ve~ v *l h » T min(l, e^-E^/kn^ (1) 

which fulfills the condition of detailed balance. As long as 
the general kinetics are concerned, characterized by the 
tendency of the system to approach thermal equilibrium, 
the specific form Q of the hopping rates should not be 
of crucial importance. 

In many cases, 25,26 including Co deposited on Pt(lll), 
adatoms on the top of a terrace or at step edges can 
exchange positions with an atom underneath in one con- 
certed move. In comparison to single atom moves consid- 
ered so far, a larger binding energy has to be overcome in 
order to effect such a simultaneous move of two particles. 
However, at the same time the transition state energy can 
be reduced because of the higher coordination in the tran- 
sition state configuration. Direct exchange processes are 
found for Co deposited on Pt(lll) over a wide temper- 
ature range 250-520 KiS 5 ^ They are especially frequent 
for low-coordinated atoms on top of terraces (with coor- 
dination 3) or at step edges (with coordinations 4 or 5). 
We therefore include direct exchange processes between 
unlike atoms in the simulation. The corresponding rate 
involves an exchange barrier U x that adds to the migra- 
tion barrier Ut, while initial and final energies are taken 
into account as before. These exchange processes, how- 
ever, are allowed only if one atom has a coordination in 
the range 3-5 and the other one in the range 8-10. Let 
us remark that the activation energy for exchange diffu- 
sion of Pt on Pt (100) is just 470mVrl which is almost 
the same as the total barrier Ut + U x when U x = 5, as 
used in most of our computations described below. 

The time evolution in our simulation is determined by 
a rejectionless continuous time MC-algorithm that gen- 
erates realizations of the master equation describing the 
growth kinetics. Each simulation step consists of: (i) in- 
crementing the time by an interval At drawn from an 
exponential distribution with mean (At) = (S/^j/) -1 ! 
(ii) executing a process i —*■ f with probability propor- 
tional to its rate Wi / , (iii) updating of rates of processes 
affected by the moving atom. 

The simulation box consists of 15 layers with 1440 lat- 
tice sites each. To generate an isolated cluster attached 
to the weakly binding substrate, we start from a seed 
with 5 atoms on the surface. For the deposition rate 
F = 3.5 ML/s used in most simulations and tempera- 
tures T > 0.5, the distance between clusters, which would 
be observed by self-organized nucleation, exceeds the box 
size. This could be effectively accounted for by increasing 
the deposition rate at the boundary of the box to match 
the extra flux of atoms deposited within the typical cap- 
ture zone of an island. However, in order to reduce the 
number of parameters, we deliberately kept the deposi- 
tion rate constant for all temperatures. 

An example of a typical cluster with N = 2000 atoms 
obtained in the KMC simulations is shown in Fig. ^ 
The top facet has (111) orientation as expected for the 
(111) substrate, and side facets are of (111) and (100) 




FIG. 1: A typical cluster configuration with 2000 atoms for 
h = 4, T = 1.2 and C/ x = 0. Pt atoms are marked in blue 
and Co atoms are marked in different colors depending on 
the number of nearest neighbor Pt atoms: if more (less) Pt 
atoms are found out of plane than in plane the atom is marked 
in green (yellow); otherwise it is marked in red. Facets with 
sixfold {111} and fourfold {100} symmetry can be clearly dis- 
tinguished. 



TABLE I: Parameters used in the simulation. All energies are 
given in units of fcsTo/1.83 = 45meV. 
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type. The surface shows strong Pt segregation. Similar 
to the experiments? 6 " the aspect ratio of lateral to height 
dimension of the cluster is about 3:1. At the temperature 
T = 1.2, the ratio D/Fa 4 ~ 2 • 10 10 , and the realization 
of the growth process implied 3.3 • 10 10 elementary pro- 
cesses. A simple MC algorithm which generates equiva- 
lent dynamics, but is not rejection free, can be estimated 
to require about N 2 D/(a 4 L 2 F) ~ 6 • 10 13 trials for ele- 
mentary moves, i.e. about 2000 times more than the al- 
gorithm used here. This ratio increases further at lower 
temperatures, since bond breaking processes will acquire 
lower rates and consequently contribute to larger time 
increments At. Subsequent equilibration of the clusters 
under zero flux during a multiple of the deposition time 
did not induce significant changes in the cluster configu- 
rations and will therefore be ignored in the following. 

To generate homogeneous films (where the width of 
the interfacial growth zone is much smaller than the av- 
erage film thickness), discussed in section iTVl we choose 
a stronger surface potential V s = —15 which guarantees 
complete wetting in the first layer in the temperature 
range considered. 



B. Structure-induced magnetic properties 

In CoPt3-alloys magnetic moments are mostly due to 
the Co atoms, with fi Co ~ \.7^ B ^2& while Pt atoms 
carry a comparatively small induced moment /i pt of 
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about 0.3/is, which has been shown to depend on the ac- 
tual atomic environment^ Hybridization of d-electrons 
between neighboring Co and Pt atoms and a strong spin- 
orbit interaction near the Pt sites lead to a magnetic 
anisotropy that tends to align the Co moment along the 
Co-Pt bond direction. In the following we adopt a bond 
picture as in earlier work on bulk systems f 2iL21i where 
the structural part of the magnetic anisotropy, Ha, is 
expressed as a sum over bonds (i, S) that connect a Co 
atom with moment fif° at site with a species a (a — 
Co,Pt,V) at a site R t + S, 

Ha = - J2J2 AC ° a ^° ■ <*)7(Im?°II<5|) 2 (2) 

<i,5) a 

By 8 we have denoted the 12 possible nearest neighbor 
bond vectors in the fee lattice. Note that this anisotropy 
term naturally entails both surface and bulk contribu- 
tions by considering vacancies as a possible neighbor 
species. 

The parameters A Coa are the anisotropy energies as- 
sociated with a Co-a bond. These are deduced from ex- 
periment and the term a — Pt gives the dominant contri- 
bution. Note that for saturated magnetization a nearest 
neighbor contribution to the anisotropic part of dipolc- 
dipole interactions has exactly the same form, so that this 
contribution is included already in the coefficients A Coa . 
The isotropic part of the exchange interactions only leads 
to a small renormalization of chemical interactions. 

Equation J5J provides a relationship between a given 
atomic structure, to be obtained from simulations, and 
the magnetic anisotropy energy. Magnetic anisotropics 
for disordered alloys, including Co-Pt alloys, have been 
obtained by microscopic calculations based on the KKR- 
CPA method^ Although less accurate in principle, the 
local expression J2J here allows us to relate effects in the 
atomic short range order to the anisotropy energy in a 
direct way. Note that 2nd order contributions of the 
anisotropy in the direction cosines as in eq. J3J) do not 
yield a magneto-crystalline anisotropy in an fee lattice 
occupied by a single species. Additional bulk contribu- 
tions are of 4th order and generally much weaker than the 
magnetic anisotropy in lower symmetry configurations to 
be considered here. 

Several observations support the bond picture under- 
lying eq. as a reasonable approximation. Magnetic 
torque and magneto-optical Kerr-effect measurements 
on Co-Pt multilayers show that the anisotropy energy 
in these systems is dominated by an interfacial contri- 
bution connected to Co-Pt bonds. Detailed measure- 
ments of multilayers yield K CoPt = 0.97 mJ/m 2 for 
(111) orientation and K GoPt = 0.59 mJ/m 2 for (100) 
orientationii^^ Considering the different angles of bonds 
to the surface and the different packings, eq. @ indeed 
can reproduce this difference with one consistent value for 
A GoPt w 250 [ieV per Co-Pt bondi^ Hence, on a semi- 
quantitative level, anisotropy energies of Co-moments in 
different chemical environments are consistent with a su- 
perposition of bond contributions and a bond energy 



^CoPt ag gj ven a b ve. 

Similarly, from measurements on Co-vacuum 
interfaces'^ and theoretical investigations of a freely 
standing Co-monolayer— we estimate A CoV ~ — 67^teV. 
The remaining parameter is A CoCo for which we retain 
only the nearest neighbor dipolar contribution as dis- 
cussed above. Using fi Co = 1.7/ab and \6\ = 2.72 A in 
CoPt 3 , this yields ,4 CoCo = 23 /xeV. 

It is now straightforward to express the part of the 
anisotropy energy caused by local chemical order, 

E s = ii^ {/-tin plane} — i?A {/a out of plane} (3) 

in terms of the numbers of Co-a bonds in plane, n Goa , 

and out of plane, n^ 00 - Using a site occupation number 
representation and the condition that occupation num- 
bers for the three species a add up to unity at each site, 
© can be brought into the form 

E s = E Goa = - (A Coa ~A GoV )P Coa 

a=Co,Pt Q=Co,Pt 

(4) 

with structural anisotropy parameters 

P Coa = ^ n ± a -^° a )- ( 5 ) 

These parameters enter as the primary structural char- 
acteristics related to the magnetic anisotropy. 

In addition one has to take into account the magnetic 
form anisotropy due to dipolar interactions. Dipolar 
sums for a given cluster with saturated magnetic moment 
M s are carried out using moments for Co and Pt as given 
before. By E^ip we denote the difference in dipolar ener- 
gies for M s parallel and perpendicular to the substrate, 
respectively. The anisotropy within the substrate plane 
turns out to be negligibly small because of the approxi- 
mate 3-fold symmetry of cluster shapes. Nearest neigh- 
bor dipole-dipole interactions were already incorporated 
in the quantity E s , as described above. The total mag- 
netic anisotropy energy of a cluster is then given by 

Etot = E s + Edip (6) 

Dipolar interactions generally favor in plane magnetiza- 
tion of thin films and clusters, Ea ip < 0. For thin homo- 
geneous films, Edip = — (^o/2)M 2 . 

As mentioned above and supported by the above esti- 
mates, the term a =Pt dominates expression J3J|. PMA 
hence requires P CoPt to be sufficiently large that E s can 
overcome the negative dipolar energy. 

III. CLUSTERS 

The two anisotropy parameters p CoPt and p CoCo de- 
fined in JSJ are shown in Fig. [21 as a function of tem- 
perature for clusters consisting of 1000 atoms. The solid 
lines connect data points for different values of the extra 
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FIG. 2: Structural anisotropy parameters P^f* (upper 
panel) and P^Tjf 30 (lower panel) and corresponding structural 
anisotropy energies. The dashed lines correspond to data ob- 
tained with a deposition rate of F = 0.35 ML/s, which is 10 
times smaller than the rate used for the full lines. 



barrier for exchange processes U x , and for a flux F = 3.5 
ML/s. In the absence of exchange processes (f7 x = oo) 
both anisotropy parameters are smaller than zero for all 
temperatures, leading to a preference of in-plane mag- 
netization. The reason is that the clusters are relatively 
flat and that surface segregation of Pt is weak because 
of kinetic suppression. As a consequence, in the surface 
of a cluster there exist more Co-Pt bonds in-plane than 
out-of-plane, while atoms in the inner part of the cluster 
give essentially no contribution to P. The flat cluster 
shape, on the other hand, would favor PMA if strong Pt 
segregation as realized at equilibrium could build up. 

Indeed, when including exchange processes, Pt segre- 
gation gets enhanced considerably and P can become 
positive. Actually, the cluster displayed in Fig. show- 
ing pronounced Pt segregation to the (111) and (100) 
facets, was computed with U x = 0. As seen from Fig. EJs,, 
for U x = 0, P CoPt changes sign near T = 0.58 and reaches 
a maximum at an optimum temperature T max ~ 0.8. 
Conversely, P CoCo < for all parameters, which leads to 
a small reduction of the sum (0} relative to its leading 
term a =Pt, see Fig. [2b- 
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FIG. 3: Concentration of Pt atoms in the outer shell of clus- 
ters for h = -Eptpt — -EcoCo = 4 in the upper panel and h = 
in the lower panel. The uppermost lines result from the mean 
field equations and ||SJ. 



The third set of data in Fig. [21 (triangles) refers to 
U x = 5. This means that the activation energy for direct 
exchange is twice the diffusion barrier. Yet exchange pro- 
cesses have an important influence as they still can render 
pCoPt p 0S iti ve The onset of a positive P CoPt and the 
temperature were P CoPt takes its maximum are shifted 
to somewhat higher values than in the case U x = 0. Fig. [21 
also contains data for the flux F = 0.35 ML/s. The re- 
duction of F by one order of magnitude apparently leads 
to a small increase of p CoPt . 

Investigating the mechanism of PMA in more detail we 
find that the sign of P CoPt is determined by two major 
factors which show opposing trends in their temperature 
dependence. These are the degree of Pt surface segrega- 
tion and the cluster shape, which we now discuss in more 
detail. 

The concentration Cpt of Pt-atoms in the outer shell 
of a 1000 atom cluster is plotted in Fig. [SJ again for 
a flux F — 3.5 ML/s and three values of t/ x . As dis- 
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cussed above, the observed degree of Pt surface segre- 
gation is generally smaller than in the equilibrium case 
because of kinetic hindrance: During growth, the time 
for attaining equilibrium through exchange and diffusion 
processes is limited due to continuous incorporation of 
newly deposited atoms. Fig. reveals that both direct 
exchange processes and an increasing temperature act to- 
wards restoring equilibrium, i.e. they facilitate Pt surface 
segregation. 

Clearly, a large Cp t in the topmost layer leads to an en- 
richment of Co-atoms in the layer underneath. Therefore, 
for oblate cluster shapes with more surface area oriented 
in the [111] direction than in directions perpendicular to 
[111], a sufficiently strong segregation of Pt to the surface 
will induce a positive sign of p CoPt . This favors PMA, 
even for a disordered structure within the cluster interior. 
Furthermore, the weakly binding van der Waals substrate 
also allows for segregation towards the substrate during 
growth, effectively doubling the available surface with fa- 
vorable orientation. 

It should be noted that segregation of the majority 
atoms (Pt) can be observed even when the parame- 
ter h — Vptpt — VcoCo controlling segregation is zero. 
This can be understood from a mean-field argument that 
counts bonds in a slab with a free (111) surface: for a ran- 
dom fee alloy structure, the exchange of a Co-atom at the 
surface with a Pt-atom in the bulk allows the system to 
lower its energy on average by 3/, so that the surface 
concentration Cpt will exceed the stoichiometric concen- 
tration, see Fig |2| By contrast, for a fully LI2 ordered 
alloy the stoichiometric composition at the surface yields 
the lowest energy. 

In order to give a more quantitative estimate of segre- 
gation, we consider a mean field model that consists of 
three completely filled (111) layers of atoms. The first 
layer is the free surface and the third one has fixed stoi- 
chiometric composition. Exchange processes are allowed 
between the outermost two layers. This model takes into 
account the almost non-existent bulk diffusion within a 
cluster, which results in an increased concentration of Co 
in the second layer impeding further exchange. The in- 
ternal energy term per atom for a given Co concentration 
in the first layer Cc ,i can be expressed in the case of 1:3 
stoichiometry as 

U = W ~h(l-4Cco,i) + h (l-C Co ,i+4Cg 0il ). (7) 
lb I 

The energy minimum at T = is attained for Cq ,i — 
1(1 — h/(2J)) clearly predicting segregation for h — 0. 
The entropy is 

1 2 

Q = Co,Pt 2=1 

Minimizing the free energy F = U — TS then yields the 
equilibrium concentration of Co and Pt in the first layer. 
The corresponding results for the Pt concentration are 
displayed in Fig. [3J Including the entropy term leads to 
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radii l\\/l± (with line as guide to the eye). Decreasing val- 
ues signify the transition from oblate to spherically shaped 
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FIG. 5: Number of atoms per layer N z represented with a 
radius R z = N z /n and averaged over 20 realizations of the 
growth process. Layer one is the bottom layer on the (111) 
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a decrease in Pt concentration in the outer shell with 
increasing temperature and restores the stoichiometric 
concentrations in the limit of high temperatures. Such 
a negative slope with temperature is reproduced in the 
simulations for h = when exchange processes are in- 
cluded, but not in the simulations for h — 4 and h = 
without exchange processes. The positive slope of simu- 
lated segregation values with temperature indicates that 
segregation differs from thermal equilibrium due to ki- 
netic limitations. This model also does not account for 
the portions on the cluster surface which are not flat. 

The second factor important for PMA is the cluster 
shape. An increasing temperature drives the shape closer 
to the equilibrium shape which is less oblate. This is seen 
from the aspect ratio l\\/l± of the gyration radii of clus- 
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FIG. 6: Scaling plot of surface E s and bulk Edip contributions 
to the magnetic anisotropy E to t = E s + Ejip as a function 
of cluster size N. The lines fitting E s and Eiip have slope 
2/3 and 1 respectively. The line through Etot is given by 
Etot = K a N 2 ' 3 — K dip N, where the prefactors K a and K to t 
are obtained from the fits to E s and _Bdi P . 

ter sizes in the direction parallel and perpendicular to the 
substrate, shown in Fig.0] Cluster shapes are parameter- 
ized in more detail in Fig. El The atomic layer z contains 
N z atoms and has an effective radius R z = y/ N z /n. The 
figure displays the connection between z and R z for a 
series of temperatures. It clearly shows the evolution of 
cluster shapes from oblate to almost a half-sphere as tem- 
perature is increased. In this way temperature effects in 
the cluster shape counteract the temperature dependent, 
segregation induced PMA. The result is the existence of 
a certain temperature window where PMA prevails, in 
agreement with experiments A£ 

The discussion so far makes it clear that PMA essen- 
tially emerges as a surface effect. This suggests that 
- nP oPt - N 2 / 3 , and accordingly P CoPt ~ iV" 1 / 3 . 
As shown by the solid line in Fig. El this behavior is well 
obeyed by the simulated data (triangles). 

Combination of our results for P Coa with (0} yields 
the structural part of the anisotropy energy. Adding 
the dipolar energy (see section |HJ we obtain the total 
anisotropy energy ©, which can be written as 

= K S N 2 / 3 - K dip N (9) 

The temperature dependence of E tot is shown in Fig. [7| 
for different t/ x and N = 10 3 atoms. An optimal temper- 
ature where E to t is maximum can clearly be identified. 
For example, for U x = 5 we have T max ~ 1. Fig. El 
contains a double-logarithmic plot of Etot as a function 
of N at T = 1, Uyc = 5. From E s and Edi P we obtain 
K s ~ 285 fieV and Edip — 5.6 ^eV. The line for E tot in 
Fig. El represents the expression © and, upon extrap- 
olation to larger iV-values, predicts the existence of an 
optimal cluster size for PMA, which is iVopt — 4 • 10 4 . 
PMA is expected to disappear when N ~ 1.3 • 10 5 . 




FIG. 7: Magnetic anisotropy energies Etot (solid lines) and 
E s (dashed lines) as a function of temperature for different 
values of the additional barrier U x for exchange processes (no 
exchange processes for (7 X = oo). E to t also includes the form 
anisotropy from dipolar interactions, while E s only contains 
the bond contributions (with inclusion of the nearest neighbor 
part of the dipolar interactions) . 



Experimental values for E to t for two different cluster 
sizes can be inferred by measurements of the blocking 
temperature using SQUID devices^ For clusters with 
N = 300 at room temperature (T = 0.56), this leads 
to the estimate E to t — 2.8 meV. Similarly, E to t = 3.4 
meV for N = 1200 at 573 K (T = 1.1). In a different 
set of experiments at T = 0.56, granular nanostructures 
(dense covering of surface, but not touching) with simi- 
lar lateral size of 3 nm were obtained with considerably 
larger anisotropy constants K ~ 13 fieV per atomi In 
our simulation model values E to t i5 12meV for N = 300 
and E tot ^ 20mcV for N = 1200 are found for T > 0.9. 
The comparison shows that the experimentally observed 
values of the anisotropy energy lie within the ranges pre- 
dicted by the model. The degree of agreement with 
experiment must be regarded as satisfactory in view of 
the uncertainties in the analysis of experiments and the 
simplifying assumptions in our model. One should also 
note the slower deposition rate in the experiments with 
F = 0.02 ML/s, compared to E = 3.5 ML/s used in most 
simulations. An effect of slower deposition rates in the 
simulations can be seen in Fig. [21 

Another effect to be mentioned is that larger spin-orbit 
couplings are known to occur for Co-atoms with low co- 
ordination on terraces or along steps on Pt-surfaccs. 40 
Such an effect in principle can further enhance PMA in 
nanoclusters because of the flat cluster shape. However, 
because of strong Pt surface segregation in the nanoclus- 
ters most of the near-surface Co-atoms are buried in the 
second layer or have high coordination. Hence we expect 
that the few Co-atoms found in the outermost layer need 
not be treated separately. 

Next we turn to analyzing the Ll2-ordering of clus- 
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FIG. 8: Dependence of the order parameter for LI2 structure 
on temperature for segregation controlling parameter h = 4 
(upper panel) and h = (lower panel). Values are averaged 
over 20 realizations for h — 4 and 10 realizations for h — 0. 



ters. Analogous to experiments, we determine the LI2 
order parameter from the magnitudes of scattered inten- 
sities around the associated three superstructure peaks 
Kj, i = 1,2,3 in reciprocal space. For that purpose, 
atoms of the cluster and the surrounding vacancies are 
represented by pseudo-spins si € {1,0,-1} at the lat- 
tice positions R; with s; = for vacancies and s; = ±1 
for A and B, respectively. The structure factor is then 
calculated from the amplitudes i*k = J^i s;e _lk Rl which 
account for both the atomic arrangement in the cluster 
and the cluster shape. The finite dimensions of the clus- 
ter lead to a significant broadening of peaks in Fourier 
space, especially in the vertical direction. In order to ac- 
count for this broadening, an integration in k-space was 
performed around each peak in form of a sphere with ra- 
dius 0.1 /a. The total scattering intensity is calculated as 
-Ki|<o.i/a l^k| 2 - As an order parameter 



adorn 



thin — Ira 



(10) 



Here Jli 2 and i ra ndom are the intensities of clusters with 
identical shape, but perfect LI2 order and random occu- 
pation by A and B atoms, respectively. 

Fig.[H]shows the temperature dependent order parame- 
ter ^> in cases of a strong surface field h = 4 (a) and h = 
(b). The inclusion of exchange processes (U x < 00) di- 
minishes \t due to reduced ordering through segregation 
effects. This can be understood from the large energy a 
Co adatom gains when it is incorporated by an exchange 
process into deeper layers of the cluster irrespective of its 
contribution to ordering. 

In the presence of exchange processes we notice from 
the figures that ordering sets in at temperatures near 
T ~ 0.8 to 1 (420 - 523 K). In experiments, the onset 
of LI2 ordering was found for temperatures near 423 K 
for 3 nm thick films consisting of a dense assembly of 
islands £ This temperature range is somewhat lower but 
close to the corresponding range in the simulations. 

Comparison with Fig. [7| shows that PMA for C7 X = 5 
starts to decrease for T > 1, which coincides with the 
onset of LI2 ordering. One should bear in mind, however, 
that PMA in clusters is a surface induced effect and its 
decrease is caused by the change in cluster shape rather 
than by the concomitant bulk ordering. 

To better visualize the statistics of growth, we present 
in Fig.[5]the layer-resolved evolution of one cluster. One 
can see pronounced incubation periods, before a nucle- 
ation event on top of the last layer takes place. Such an 
event is accompanied by a rapid initial growth and com- 
pletion of the newly formed layer, followed by a period 
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FIG. 9: Atoms per layer versus total number of atoms during 
the growth of one cluster. The curve on top shows the number 
of atoms in the first layer above the substrate and curves 
below represent successive layers on top. For the 4th to 8th 
layer nucleation events are visible as a strong initial increase 
in the number of atoms. 
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FIG. 10: Anisotropy parameter Pbuik for films as a function 
of temperature for different values of the exchange energy 
barrier (7 X and surface field h. The error bars are calculated 
from the statistical error after averaging over 5 realizations of 
the growth process and the fit uncertainty (P as a function 
of 1/N, see text) with error propagation. 

of essentially lateral growth of the whole cluster with- 
out change in height. The slope in the atom numbers of 
the top layer immediately after nucleation is much higher 
than one. This means that most of the atoms incorpo- 
rated in a new layer during its completion arrive through 
mass transport along the side facets, leading to a tran- 
sient lateral shrinkage of all the layers below the growing 
top layer. This is seen in the figure by the small dips in 
all curves occurring simultaneously during the short time 
intervals when the top layer is filled. 

The nucleation events on top can be described by a 
theory originally devised for second layer nucleation^ 1 " 
This theory describes nucleation in a confined geometry 
with an influx of particles through deposition and loss 
of particles over a step edge. An essential parameter is 
the Schwoebel barrier that atoms have to surmount when 
crossing the step edge. In our simulations, this extra en- 
ergy barrier effectively is about AE S — Vo because the 
intermediate state for a transition from the top facet to 
a side facet has one bond less than adatoms with coor- 
dination three on the top surface. With the parameters 
typically used in the simulations, the fluctuation domi- 
nated regime III discussed in Ref. ^3 1S the relevant one 
for critical nuclei of size i = 1,2. However, with the ex- 
tra influx of atoms from the side facets, it is likely that 
the mean number of adatoms on top of the facet satisfies 
n > i, so that mean-field theory may become justified for 
top layer nucleation^i 



IV. FILMS 

The same model is now applied to continuous films. 
We simulate the growth of films starting with an empty 



FIG. 11: Long range order parameter ^ for films as a function 
of temperature for different values of the exchange energy 
barrier [7 X and surface field h. Values are averaged over 5 
realizations of the growth process. 



simulation box with 224 atoms per layer and simulate 
growth up to deposition of 29 layers. The surface binding 
energy of V s — 3Vq = — 15 induces film instead of cluster 
growth. As before, the surface binding is independent 
of the atom type so that segregation to the substrate 
similar to a free surface can be observed. Consequently, 
true bulk properties are only found a few layers above the 
substrate and below the top surface. In order to separate 
the bulk contribution from surface effects, we considered 
each realization of the growth process in the range of 
2000 to 6500 atoms and made a linear fit of P vs. 1/N. 
The anisotropy parameter in the bulk -Pbuik is then found 
by extrapolating N — ► oo, 43 Fig. I1UI shows the temper- 
ature dependence of the anisotropy parameter Pbuik for 
different values of U x and h. For most sets of parameters, 
-Pbuik < 0, i.e. Co-Pt bonds align preferentially in the film 
plane. However, for certain combinations of parameters, 
the simulations yield Pbuik > 0, supporting PMA. 

Clearly, even if Pbuik > it is too small to account 
for the measured PMA in thick CoPt3 films at elevated 
temperatures. As shown recently, those measurements 
can be explained by a model with interatomic interaction 
potentials that depend on the coordination of the atoms 
involvediii One essential feature of this model is that 
an effective segregation parameter, related to h in our 
model, changes sign for intermediate coordinations of Co- 
atoms and thus favors Co segregation towards step edges 
in the surface. This sign change of h, however, should 
not significantly affect our conclusions. There we focused 
on the temperature range T < 1, where Co-segregation 
according to Ref. ceases to be effective. 

The degree of long-range order in films is shown in 
Fig. llll at temperatures below the transition temperature 
To ~ 1.83, where bulk kinetics are suppressed. With 
higher T, the order parameter ^ becomes larger. Similar 
to our findings for clusters, the order parameter 'J in the 
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case h = 4 is suppressed relative to h = 1, as clearly seen 
from the data at T — 0.8. Surface segregation therefore 
appears to impede the development of long range order. 

V. CONCLUSIONS 

In order to study growth kinetics of binary alloys we 
developed a lattice model based on nearest neighbor 
bonds with bond energies chosen to match equilibrium 
properties of CoPt3. The kinetic parameters were ob- 
tained from diffusion experiments^ DFT calculations^ 
and in addition from experiments observing interatomic 
exchange processes in the Co/Pt systemi2S4i 

Experiments on CoPt3 nanoclusters on a weakly bind- 
ing substrate^ have revealed PMA in a temperature win- 
dow that is well reproduced by our simulations. This 
temperature range is bounded towards low temperatures 
by frozen-in surface kinetics. The disappearance of PMA 
at higher temperatures is explained by the interplay of 
Pt surface segregation facilitated by direct exchange pro- 
cesses, and a transition from oblate to spherical cluster 
shapes. Our analysis suggests that the transition is not 
caused by LI2 ordering. Yet, in our simulations the on- 
set of LI2 ordering is detected in the same temperature 
range where PMA disappears, in qualitative agreement 
with the measurements. It should be noted that up to 
T ~ 1 long range order is induced solely by surface pro- 
cesses, as the bulk kinetics remain still frozen. 

The structural anisotropy responsible for the magnetic 
anisotropy is characterized in our model primarily by a 
difference in the numbers of Co-Pt bonds out-of-plane 
and in-plane. Within a bond picture each Co-Pt bond 
provides a local contribution to the magnetic anisotropy, 
which tends to align the Co-moment parallel to the bond. 
An associated magnetic bond energy, deduced from ex- 
periments on Co-Pt multilayers, can reasonably repro- 
duce the magnitude of PMA measured for nanoclusters 
in the appropriate temperature range. PMA is predicted 
to be a surface effect, a feature which could be tested 
experimentally. 

Dipolar interactions tend to turn the easy axis of mag- 
netization into the plane. The two competing effects, 



bond anisotropy (surface term) and dipolar interactions 
(volume term), lead to an optimum cluster size where 
the anisotropy energy of a cluster is largest, and a sec- 
ond characteristic size where PMA disappears. 

The importance of cluster shape effects revealed by our 
computations is also corroborated by experiment, where 
a rotation of the easy magnetization axis into the film 
plane and an associated change in the aspect ratio from 
oblate to prolate was found in separate measurements on 
cluster assemblies at higher coverages £ In our simula- 
tions, we observed clusters to become more spherical in 
shape and to develop an excess of in-plane Co-Pt bonds 
with increasing temperature. 

Application of the same methodology to growing con- 
tinuous films yields a bulk structural anisotropy favoring 
in-plane magnetization for most of the parameters stud- 
ied. However, it can yield a positive bulk contribution 
to p CoPt in a certain temperature range provided direct 
exchange processes are sufficiently fast. The associated 
PMA, however, is quite small and disappears for T > 1. 
At these higher temperatures a different mechanism has 
been proposed recently to describe PMA in thick films in 
terms of in-plane Co clustering^ 

Finally, we like to draw attention to the question of 
growth of magnetic clusters or films in strong external 
magnetic fields that saturates the magnetization in the 
growth direction. The magnetic anisotropy energy should 
then induce an anisotropy in the probabilities for atomic 
hopping such that out-of-plane Co-Pt bonding and PMA 
become favored. In fact, simulation results for growth 
in a magnetic field suggest that this effect might be- 
come detectable in the Co-Pt system. In particular, for 
Llo ordered alloys estimates based on Landau theory 
suggest that this additional field-induced PMA becomes 
stronger 4& 
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